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A new scheme for constructing approximate effective electron potentials 
within density-functional theory is proposed. The scheme consists of calcu- 
lating the effective potential for a series of reference systems, and then using 
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these potentials to construct the potential of a general system. To make 
■ contact to the reference system the neutral-sphere radius of each atom is 

used. The scheme can simplify calculations with partial wave methods in the 
q\ \ atomic-sphere or muffin-tin approximation, since potential parameters can be 

precalculated and then for a general system obtained through simple interpo- 
lation formulas. We have applied the scheme to construct electron potentials 
of phonons, surfaces, and different crystal structures of silicon and aluminum 
atoms, and found excellent agreement with the self-consistent effective poten- 
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Q , tial. By using an approximate total electron density obtained from a superpo- 

sition of atom-based densities, the energy zero of the corresponding effective 
potential can be found and the energy shifts in the mean potential between in- 
^ ■ equivalent atoms can therefore be directly estimated. This approach is shown 

to work well for surfaces and phonons of silicon. 
One route that seems promising in order to construct computationally efficient 'ab ini- 
tio' schemes for calculating total energies and forces of solids is to exploit the variational 
properties of density-functional theory [[I]]. We have shown earlier how the total electron 
density can be decomposed into a superposition of transferable atom-based densities for 
metals and semiconductors [@J3||. When such densities are used to generate an input density 
for the Harris functional excellent total energies are obtained for surfaces, phonons 

and structural differences |||3| due to the fact that the Harris functional is stationary in the 
density. In this way the self-consistency loop is avoided. It is the purpose of the present 
report to show how the the variational nature of density-functional theory can be exploited 
even further by working with both approximate densities and potentials simultaneously. 

The Hohenberg-Kohn density functional can be generalized to a functional E[n, v] which 
depends on both the density n and the potential v [|5],|7J and which is stationary with respect 
to independent variations of the density and the potential. The general functional can be 
written 
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E[n,v]=^2e a [v]- n(r)v(r)dr + E el [n] + E xc [n], (1) 

where e a denotes the eigenvalues generated by the potential v, and where E e i[n] and -ErcH 
is the electrostatic and exchange-correlation energy functionals, respectively. If the potential 
is restricted to be a functional of the density, the Hohenberg-Kohn functional or the Harris 
functional appear as special cases J7|. The stationary property of the general functional 
with respect to variations in the potential can be utilized to construct efficient schemes 
for evaluation of total energies and Hellman-Feynman forces 0. In the following we shall 
describe one such scheme which has its root in the effective- medium theory 0. The scheme 
applies to situations in which the kinetic energy can be calculated within the muffin-tin 
or atomic-sphere approximation (ASA) with spherically-symmetric potentials within the 
atomic spheres. The idea is to use self-consistently calculated potentials from a series of 
reference systems which we choose here to be a bulk crystal with varying lattice constant. 
For a given atom in a general system the potential within the atomic sphere around the 
atom is then approximated by the potential in the reference system with an appropriate 
lattice constant. The lattice constant of the associated reference system is determined by 
the requirement that a neutral sphere around the atom should have the same radius in the 
system under study and in the reference system. 

The neutral-sphere is a sphere containing 3 (4) electrons in the case of aluminum (sili- 
con) in the pseudo-potential scheme. If the approach described here is combined with the 
density construction of Ref. || where the total electron density, n(r), is approximated by a 
superposition of atom-based densities, An aiom , positioned at each atomic site, Rj, 

n(r) =^An atom (|r-R i |), (2) 

i 

then the NS radius can be obtained directly, or simple interpolation formulas can be made 
from which the NS radius can be obtained with high accuracy |]J| . 

We have used this scheme to calculate the effective potential of silicon and aluminum 
atoms in different configurations. For silicon we use the diamond structure and for aluminum 
the fee structure as reference systems. The lattice parameter is regarded as a parameter 
which can be varied in order to find a good approximation for the potential. To calculate 
the effective potentials we use a self-consistent plane-wave pseudo-potential program, with 
a 12 Rydberg cutoff for the plane-wave basis set. With this cutoff, we get the lowest energy 
configuration of silicon to be the diamond lattice with lattice constant 10.17 ao, and for 
aluminum the fee lattice with lattice constant 7.48 ao- 

As test systems for silicon we consider the diamond longitudinal phonon at the X point 
(denoted LAO(X) and frozen at the displacement 0.02 in units of the lattice constant), 
the diamond (100) surface, and the fee structure with a lattice constant of 7.18 ao- And 
similarly for aluminum, we use the fee longitudinal phonon at the X point (denoted L(X) and 
displacement 0.02), the fee (100) surface, and the diamond structure with a lattice constant 
of 11.05 ao- These six structures cover the two elements in quite different surroundings, 
and if the potential construction works for these situations a high degree of transferability 
is guaranteed. 

As the first test system, we consider silicon in the fee structure. As noted in the in- 
troduction the way we make contact between the test system and the reference system is 



2 



through the neutral-sphere (NS) radius. In the fee structure the NS radius is almost equal 
to the Wigner-Seitz (WS) radius, while the NS radius is substantially smaller than the WS 
radius in the diamond structure. This difference is due to the large regions in the diamond 
lattice which contain almost no charge and therefore are not included in the neutral sphere. 
In Fig. 1 we show the local part of the self-consistent effective potential of the fee test system 
compared to that of the reference system with the same NS and WS radius as the fee test 
system, respectively. Clearly the potential of the reference system chosen according to the 
NS criterion gives by far the best approximation to the fee potential. 

In order to quantify the difference in the potentials, we introduce the r.m.s. error a, 
defined by 



where s w is the WS radius of the test system, and the two potentials are aligned such 
that they have the same average within the sphere. To get an estimate of the error in the 
potentials due to the finite plane- wave basis set, we have compared one of the reference 
potentials of silicon with that of an 18 Rydberg calculation. We find a r.m.s. error of 
a = 0.06 eV, so this is the level of accuracy we ideally can obtain. 

In Fig. 2 we show the r.m.s error between the potential of the fee test system and the 
diamond reference system as a function of the WS radius (i.e. as a function of lattice 
constant) of the diamond lattice. We see that optimally the reference system should be 
chosen with a WS radius of 3.18 (Xq. We note that this is not at all close to the WS radius 
in the fee test system (s w = 2.806 ao). In the ASA the reference diamond structure is often 
embedded in a bec structure with twice the number of spheres where only half of them 
contain an atom, and the others are empty. For this construction the sphere radius is a 
factor 2 1//3 smaller than the WS radius, and the reference system where this smaller sphere 
equals the fee WS radius is given by s w = 3.53 ao, which still is far from the the optimal 
reference system. However the reference diamond lattice with the same NS radius as the fee 
test system, has a WS radius of 3.178 ao, and is therefore almost exact the optimal reference 
system. This is not just a coincidence, since for all test systems investigated we have found 
this to be the case. 

In Table 1 we show the minimal r.m.s. error between the potential of the six test systems 
and the reference system, compared to the error when the reference system is chosen to have 
the same NS radius or the same WS radius as the test system. It is evident from the table 
that using the reference potential chosen according to the NS criterion is almost optimal, 
while the WS criterion is far from optimal. It should be noted that if the integral in Eq. (2) 
were done within the neutral sphere instead of within the WS sphere, the NS criterion would 
give even smaller errors. 

In order to estimate how much the errors induced by the approximations for the potential 
will affect the total energy, we show in Table 2 the value of LMTO potential parameters 
H for the three potentials of Fig. 1. The potential parameters are calculated by solving 
the radial Schrodinger equation at a fixed energy for the s and p angular component (e u = 
{13.22, 18.96}), with the energy chosen at the center of gravity of the occupied part of the 
fee band. The accuracy of the potential parameters directly reflects the accuracy of the 
corresponding one-electron bands, and thereby the one-electron band energy. As seen from 
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Table 2 the potential parameters obtained with the NS criterion are in excellent agreement 
with the self-consistent parameters, while the potential parameters obtained with the WS 
criterion are more than 10 percent off. 

Until now we have aligned the average potentials within the WS sphere. However, for 
a surface calculation the mean potential at the surface is shifted relative to the bulk and 
we need to describe this shift in order to get a good estimate of the overall potential. At 
first sight the problem seems hard to overcome since the mean potential is arbitrary for 
a bulk calculation due to the divergence of the electrostatic potential ||. However it is 
possible to circumvent this problem if we use the potential constructed from a superposition 
of atom-based densities. We will name this potential the Harris potential since this is the 
effective potential used in the Harris functional, when the input density is obtained from 
a superposition of atom-based densities. We know from previous studies that when this 
potential is used as an input to the Harris functional, excellent total energies are obtained 
for phonons, surfaces, and different crystal structures 

The electrostatic part of the Harris potential V e f arrts can in a natural way be divided 
into a sum over atom-based electrostatic potentials v e i(r) each given by the sum of one ionic 
potential Vi on and the Hartree potential derived from one atom-based density An atom 

K i (r)=E^(|r-R,|)=EK n .(|r-R ? |)+ / An ^^~ R ^ rf r '). (4) 

i i 'II 

With this construction the 1/r divergence of the electrostatic potential is avoided, since the 
atom-based density decreases exponentially, and thereby fixes the vacuum level. Note that 
a consequence of this is that within this approximation all surfaces of a solid will have the 
same work function. In contrast to the total energy the work function is not variational 
in the density, and we cannot expect the density ansatz Eq. (3) to produce an accurate 
estimate. 

Having established a common energy zero for all Harris potentials we can now proceed 
to determine the energy shifts of the mean potential for atoms in different environments. In 
Fig. 3 the solid curve indicates the mean Harris potential for silicon in the reference diamond 
structure as a function of the NS radius. Also shown are the actual shifts of both the Harris 
and the self-consistent potential at the three principal surfaces and the potential shifts of 
the two inequivalent atoms in the LAO(X) phonon. These potential shifts are marked in 
the figure at the calculated NS radii of the atoms. We see that the potential shift in the 
reference system compares surprisingly well to the potential shifts in the test systems. On 
the average, the shift of the Harris potential is about 4 percent higher than for the self- 
consistent potential, and the shift of the reference system is about 8 percent higher than for 
the Harris potential. 

In Table 3 we show the r.m.s. error between the Harris potential and the self-consistent 
potential for a silicon atom at the diamond (100) surface. This is compared to the r.m.s 
error between the self-consistent surface potential and both the self-consistent and Harris 
potential of the reference system when the reference system is chosen according to the NS 
criterion. As seen from the table the r.m.s errors are similar for the three potentials, and 
since the total energies obtained using the Harris potential are excellent ||[|, all three 
reference potentials are accurate enough to be used to calculate total energies. 

In summary we have presented a scheme for obtaining transferable ASA potentials. 
The scheme was applied to six test systems consisting of silicon or aluminum atoms. The 
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obtained potentials were in very good agreement with the actual self-consistent potentials, 
reproducing both the radial variations and the shifts in the mean potential, and when used 
as input to the LMTO method accurate potential parameters were obtained. We expect 
the method to be valuable for constructing new approximate total energy schemes, and it is 
presently used in a new formulation of an Effective- Medium Tight-Binding model for Silicon 
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FIGURES 



FIG. 1. The figure shows the self-consistent effective potential (i.e. r 2 v(r)) of silicon in the fee 
structure (solid line) and that of the diamond reference system with the same NS radius (dotted 
line) and WS radius (dashed line), as a function of the radial distance The potentials have been 
aligned such that the mean potential within the WS sphere is zero. The two vertical lines show 
the NS and WS radius of the fee system, respectively. 

FIG. 2. The figure shows the r.m.s error (Eq. (2)) between the effective potential of silicon 
in the fee structure and the potential in the diamond reference system as a function of the WS 
radius of the reference system. The three crosses show the error when the reference system has the 
same WS radius, NS radius and bec WS radius (diamond with empty spheres) as the fee system, 
respectively. 

FIG. 3. The solid curve shows the mean Harris potential, within the NS, of the diamond 
reference system as a function of the NS radius. The two first crosses show the shift in the mean 
Harris potential for the two inequivalent atoms in the LAO(X) phonon, at their NS radii. The 
last three crosses show the shift in the mean potential at the diamond (111), (110), and (100) 
surface, respectively. The circles show the shifts for the corresponding self-consistent potentials. 
The dotted line shows the mean potential in the equilibrium diamond lattice. 
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TABLES 

TABLE I. Calculated LMTO potential parameters for the three pseudopotentials of Fig. 1. For 
each angular component the appropriate nonlocal contribution to the potential has been added. 



potential 


C p - C s [eV] 


A s [eV] 


A p [eV] 


SC fee 


10.968 


1.566 


1.358 


ref. NS 


10.961 


1.565 


1.355 


ref. WS 


13.028 


1.786 


1.432 



TABLE II. The first row shows the minimal r.m.s. error between the potentials of the six 
test systems and the reference system. The second and third row show the r.m.s. error when the 
reference system is that with the same NS and WS radius as the test system, respectively. The 
first three columns show the errors for the three silicon test systems; the fee structure, the (100) 
diamond surface, the diamond longitudinal phonon at the X point (LAO (X)) with a displacement of 
0.02 in units of the lattice constant. The last three columns show the errors for the three aluminum 
test systems; The diamond structure, the fee (100) surface, and the fee longitudinal phonon at the 
X point (L(X)) with displacement 0.02 in units of the lattice constant. 





FCC 


Silicon 
(100) 


LAO(X) 


Diamond 


Aluminum 
(100) 


L(X) 


g.min 


(eV) 


0.15 


0.08 


0.06 


0.11 


0.05 


0.028 


a NS 


(eV) 


0.15 


0.20 


0.08 


0.24 


0.10 


0.034 


a ws 


(eV) 


1.49 


1.22 


0.34 


0.89 


0.48 


0.042 



TABLE III. The r.m.s. error of potential differences (Eq. (2)) for a silicon atom at the (100) 
surface. The differences are between the self-consistent surface potential V sc , the Harris potential 
yHams £ or sur f ace w ith the density construction Eq. (3), and the analogous potentials in the 
diamond reference system chosen according to the neutral sphere (NS) criterion. 





ySC y Harris 


ySC _ ySC 
ref 


t/SC ^/Harris 
V ~ V re.f 


a NS (eV) 


0.24 


0.20 


0.30 
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1.8 




0.0 I 1 1 1 1 1 1 1 1 1 

2.7 2.9 3.1 3.3 3.5 3.7 

Wigner-Seitz radius (a ) 



-22.5 




neutral-sphere radius (a ) 



